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1 .   INTRODUCTION 

The  study  of  atmospheric  diffusion  can  be  simply  defined  as 
the  investigation  of  the  diluting  effects  of  the  atmosphere  on  a 
released  contaminant.   The  atmosphere  is  most  often  in  a 
turbulent  state  implying  that  its  motions  can  be  predicted  only 
in  a  stochastic  sense.   Since  virtually  all  gaseous  pollutants 
and  most  airborne  particulates  very  quickly  become  well  mixed  on 
the  molecular  level  in  the  troposphere,  diffusion  of  those 
materials  can  be  described  in  a  similar  probabal ist ic  manner. 
Introducing  statistics  automatically  requires  the  user  to 
consider  the  variations  of  those  statistics  in  time  and  space. 
As  an  example  of  the  later,  a  predicted  crosswind  concentration 
distribution  represents  only  a  first  order  approximation  to  the 
actual  distribution,  which  will  often  show  significant  variations 
from  that  prediction.   The  study  of  those  fluctuations  is 
presently  only  beginning  (see  Sawford,  1985). 

This  study  focuses  on  the  time  variations  of  diffusion 
statistics,  specifically,  the  standard  deviation  of  the  crosswind 
concentration  distribution  (sigma-y).   Skupniewicz  and  Schacher 
(1984b)  have  shown  that  overwater  releases  of  material  will 
diffuse  in  the  crosswind  (y)  direction  with  two  scales  of  motion. 
On  one  scale,  material  will  disperse  about  the  centerline  of  the 
plume  or  cloud  due  to  turbulence  of  length  scales  close  to  the 
plume  or  cloud  size.   This  scale  of  dispersion  is  referred  to  as 
relative  dispersion,  implying  that  the  statistics  can  be  observed 


by  the  motions  of  parcels  relative  to  each  other  in  a  moving 
frame  of  reference.   On  the  second  and  much  larger  scale,  the 
plume  or  puff  will  disperse  relative  to  a  fixed  axis,  usually- 
chosen  to  be  the  mean  wind  direction.   Turbulence  acts  on  the 
instantaneous  plume  or  puff  as  if  it  were  a  single  entity,  and  is 
referred  to  as  single-particle  diffusion.   A  more  common  and 
understandable  synonym  is  meander. 

Skupniewicz  and  Schacher  (1984a)  implicitly  considered  the 
combined  diffusion  properties  of  both  components  when  hourly 
averaged  plume  parameterizations  were  derived  from  tracer 
experiments.   In  this  report,  the  two  components  are  separated  in 
order  to  supply  more  information  on  the  nature  of  the  diffusion 
processes,  and  resultantly  more  knowledge  of  potential  hazard 
from  a  contaminant  release. 


2.   EXPANDED  DATA  BASE 

Most  descriptions  of  plume  parameters  rely  on  empirical 
formulae  and/or  major  assumptions  about  the  physical  behavior  of 
the  atmosphere.   Verification  of  these  parameterizations 
therefore  should  utilize  data  gathered  under  a  wide  variety  of 
atmospheric  conditions  and  geographical  locations.   The 
Environmental  Physics  Group  (EPG)  at  the  Naval  Postgraduate 
School  (NPS)  has  made  a  significant  effort  to  collect  a  diverse, 
overwater  data  base  over  a  several  year  period  for  this  purpose. 
These  data  are  reported  in  Skupniewicz  and  Schacher  (1984a)  and 
will  be  hereafter  referred  to  as  the  NPS  data.   The  NPS  data 
primarily  consists  of  concentration  profiles  collected  from  an 
aircraft  platform  at  various  distances  from  a  continuous  surface 
release  of  SF6  gas.   Part  of  the  work  reported  here  was  to 
integrate  a  new  data  set  collected  by  the  German  Military 
Geophysical  Office  (GMGO)  (see  Groll  et  al,  1984)  into  our  data 
base.   These  data  were  obtained  in  the  North  Sea  via  continuous 
SF6  releases  from  a  ship  and  subsequent  downwind  plume  transects 
with  a  second  ship  sampling  gas  concentrations.   Table  1  lists 
information  on  the  three  experiments  made  available  to  EPG. 


Experiment  Data    No.  of  Sampling  Boats    No.  of  PI 
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Table  1.   GMGO  overwater  tracer  ex^peMmehtsQVa'i  lable-'S-t  EPG-NPS. 

Each  transect  produces  an  "instantaneous"  sigma-y  value  used  in 
this  analysis. 


Task  I  specifically  involved  several  steps  needed  to  produce 
plume  parameters  from  the  GMGO  data: 

1)  transferring  original  9-track  data  files  to  mass  storage 
files  in  IBM  3033  readable  format, 

2)  combining  half-hour  average  meteorological  files  with 
nearly  coincident  tracer  profiles,  and 

3)  calculating  the  second  moment  of  the  concentration 
distribution  (sigma-y)  for  each  profile  and  the 
downwind  distance  from  navigational  information. 

Table  2  lists  the  resulting  data  set  contents  and  format. 


DEFINITION 

day,  month,  year 

hour 

relative  humidity 

wind  direction 

dry  bulb  at  17m 

wind  speed  (MET) 

wet  bulb  at  17m 

water  temperature 

dry  bulb  at  3m 

wet  bulb  at  3m 

gradient 

speed  of  boat 

source  strength 

sampling  height 

wind  speed  (SF6) 

distance  from  source 

bearing  toward  source 

heading 

no.  of  points  used  for  CWCI 

cross  wind  concentration 

integration 
mean  mass  position 
standard  deviation 

Table  2.   Contents  of  new  GMGO  overwater  tracer  data  set  added  to 
the  EPG  overwater  data  base.   Format  is  in  FORTRAN  code.   Public 
mass  storage  data  set  name  is  "MSS. F3896 . DATA. NEWSM" . 


VARIABLE 

DATE 

TIME 

RH 

WD 

DB17 

WSM 

WB17 

ET 

DB3 

WB3 

GRAD 

SB 

SS 

SH 

wss 

DS 

BTS 

HD 

#  pt 

CWCI 

MEAN 

SD 

UNITS 

FORMAT 

no  units 

312 

Z 

IH 

% 

F5.1 

deg 

m 

C 

F5.1 

m/s 

F5.2 

C 

F5.1 

C 

F5.1 

C 

F5.1 

C 

F5.1 

C/1000m 

14 

m/s 

F5.2 

m^/hr 

F2.0 

m 

12 

m/s 

F5.2 

m 

F8.1 

deg 

F5.0 

deg 

F2.0 

no  units 

12 

ppb-m 

F8.0 

m 

F6.  1 

m 

F6.  1 

3.   RELATIVE  DIFFUSION  PARAMETERIZATION 

The  basic  approach  used  in  parameterizing  the  relative 
diffusion  is  to  group  observations  according  to  Pasqui 11-Gif f ord 
equivalent  stability  classes  (see  Turner,  1967),  and  then  regress 
sigma-y  versus  range  (from  the  release  point)  for  each  class. 
Sigma-y  was  defined  in  the  last  section  as  the  second  moment  of 
the  concentration  distribution  for  a  given  profile.   The 
customary  assumption  that  concentration  in  the  y  direction  is 
independent  of  concentration  in  either  the  x  or  z  direction  leads 
to  the  conclusion  that  these  instantaneous  "snapshots"  of  the  SF6 
plume  are  identical  to  the  crosswind  dimensions  of  a  hypothetical 
"puff"  or  "burst"  release  under  the  same  atmospheric  conditions. 

The  methodology  for  determining  stability  class  is  described 
in  detail  in  Skupniewicz  and  Schacher  (1984a).   Briefly,  the 
scheme  uses  windspeed,  air-sea  temperature  differences  and 
relative  humidity  to  produce  stability  classes  equivalent  to  the 
well-known  Pasquill-Gif f ord  diffusion  categories.  [Whether  a  true 
equivalence  exists  is  a  matter  of  controversy  and  is  discussed  in 
Skupniewicz  and  Schacher  (1984a).] 

Under  this  scheme  classes  are  assigned  a  letter  designating 
atmospheric  stability,  with  "B"  representing  the  most  unstable 
situations,  "C"  representing  moderately  unstable,  "D"  representing 
neutral,  and  "E"  representing  moderately  stable  conditions. 

After  dividing  the  data  according  to  stability,  the  subsets 
were  examined  for  dependence  upon  wind  speed,  time  of  day,  and 
day  to  day  trends.   No  obvious  dependences  were  found.   Figures 


1-3  show  sigma-y  vs.  downwind  distance,  with  the  integer 
truncation  of  the  wind  speed  used  as  the  point  designator. 
In  as  much  as  windspeed  gives  a  first  order  approximation  to 
sea-state,  increased  surface  stress  within  a  given  stability 
class  has  a  negligible  effect  on  diffusion.   This  result  is 
particularly  surprising  for  the  neutral  case  (class  D) ,  where  a 
small  air-sea  temperature  difference  will  always  produce  neutral 
stability  under  this  scheme  regardless  of  windspeed.   Intuitively, 
one  would  expect  that  higher  windspeed  would  enhance  the  plume 
spread.   Because  stable  or  unstable  profiles  tend  towards 
neutrality  as  windspeed  increases  in  this  stability  scheme, 
windspeed  dependence  is  not  expected  in  non-neutral  categories. 
It  is  also  noteworthy  that  stable  situations  are  found 
exclusively  in  light  wind  situations  in  the  North  Sea  experiments. 
These  situations  were  not  found  in  the  California  experiments. 
Based  on  the  above  data,  the  parameterization  chosen  for 
relative  diffusion  is  a  simple  linear  relationship.   The 
following  factors  influenced  this  choice.   Instantaneous  releases 
into  a  turbulent  atmosphere  will  theoretically  experience  the 
following  four  different  growth  regimes: 

a)  the  initial  stage 

b)  the  inertial  stage 

c)  the  central  stage 

d)  the  final  stage 
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Figure  1.   Sigma-y  relative  vs.  downwind  distance.   Symbol  of 
data  point  is  integer  truncation  of  windspeed  in  m/s.   Only 
data   3200  m  from  source  is  plotted. 

a)  Pasquill-Gif f ord  stability  class  B  (moderately  unstable), 
GMGO  data.   23  data  points  are  hidden. 

b)  Pasquill-Gif ford  stability  class  C  (slightly  unstable), 
GMGO  data.   15  data  points  hidden,  1  data  point  missing. 
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Figure  2.   Same  as  figure  1,  except  a)   Pasquill-Gif f ord  stability 
class  D  (neutral),  GMGO  data.   194  data  points  hidden,  114  data 
points  missing.   b)   class  E  (slightly  stable),  GMGO  data.   8  data 
points  hidden,  153  data  points  missing. 
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Figure  3.   Same  as  figure  1,  except  a)   Pasquill-Gif f ord  stability 
class  D  (neutral),  NPS  data.   13  data  points  hidden,  145  data  points 
missing.   b)   class  E  (slightly  stable),  NPS  data.   1  data  point 
hidden,  45  data  points  missing. 


10 


The  initial  stage  applies  to  distances  of  only  10-100m.   The 
inertial  stage  will  be  relevent  from  100m  to  less  than  1km.   The 
central  stage  will  apply  from  1km  to  10km  and  the  final  stage 
will  apply  thereafter.   These  bounds  are  highly  variable. 

The  initial  stage  of  growth  is  valid  when  the  instantaneous 
puff  "remembers"  its  initial  size.   In  this  regime  cloud  size 
will  increase  linearly  with  distance  or  travel  time  because  each 
particle  in  the  cloud  will  move  with  its  initial  velocity  and  the 
particles  have  not  had  time  enough  to  become  correlated  with  each 
other . 

The  initial  stage  is  followed  by  growth  dominated  by 
turbulence  in  the  inertial  subrange.   In  the  inertial  region, 
plume  growth  can  be  shown  to  follow 


2      ^3 
o   =  aet 
yr 


(1) 


where  0yr  is  the  relative  sigma-y, 

a   is  a  constant, 

e   is  the  turbulent  kinetic  energy  dissipation  rate,  and 

t   is  travel  time  of  the  cloud. 

This  inertial  stage  of  growth  is  valid  only  where  the  width  of 
the  cloud  is  small  compared  to  the  average  height  of  the  cloud; 
surface-released  clouds  seldom  meet  this  criterium. 
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Smith  and  Hay  (1961)  laid  the  foundation  for  considering 
cloud  spread  as  a  function  of  turbulent  energy  in  a  "sliding" 
spectral  window.   They  found  both  experimentally  and 
theoretically  that  cloud  spread  will  approximately  follow  a 
linear  growth  relationship 


cr  =0.22  ix  , 


where  i  is  the  total  turbulence  intensity,  and 


1  = 


,,   2  2  2.v  y2 

L<u_+v_+w_>)   2 


where    x  is  downwind  distance, 

u,v,w  are  fluctuating  wind  components,  and 
U  is  mean  windspeed. 


(2) 


(3) 


Mikkelsen  and  Eckman  (1984)  have  recently  given  support  to  these 
results  from  overland  experiments  for  surface  releases  in  short 
to  medium  ranges. 

At  long  distances  from  the  source,  correlation  between 
hypothetical  particles  within  the  cloud  approaches  zero,  and  the 
spread  of  the  cloud  will  behave  like  the  asymptotic  single- 
particle  solution, 

i/ 

(4) 


0    oc   t  V* 


yr 
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Large  scale  motions,  however,  are  often  organized  into  coherent 
vortex  structures  (i.e.  cyclones).   As  a  result,  this  limit  is 
rarely  reached,  and  growth  continues  along  a  more  or  less  linear 
asymptote . 

All  of  these  above  arguments  and  the  rather  linear  shape  of 
the  data  scatter  lead  us  to  a  linear  paramterization.   The 
results  are  shown  in  figures  4-6  and  summarized  in  table  3.   Note 
that  there  is  very  little  unstable  data  at  ranges  greater  than 
3000m.   This  is  primarily  due  to  the  fact  that,  under  these 
conditions,  diffusion  is  enhanced  and  tracer  concentrations  are 
greatly  reduced  far  from  the  source.   This  makes  the  plume 
difficult  to  locate,  and  may  bias  measurements  towards  narrow, 
concentrated  profiles.   The  most  striking  feature  is  that  there 
is  no  clear  distinction  between  P-G  classes  B  and  C  ( GMGO  data), 
or  between  classes  D  and  E  (both  data  sets).   Those  pairs  are  so 
closely  matched  that  we  suggest  only  a  two-class  parameterization 

(also  shown  in  table  3) ;  one  for  all  unstable  conditions  and  a 
second  for  neutral  to  stable  stratification.   While  turbulence 
intensity  is  not  explicitly  used  in  the  paramterization,  this 
two-class  system  is  supported  by  the  calculated  wind  variances 
from  the  NPS  data  shown  in  table  4.   These  wind  data  include  some 
measurements  from  periods  not  coincident  with  tracer  releases, 
and  do  not  include  some  periods  when  tracers  were  released. 
Nonetheless,  the  largest  change  in  turbulence  (when  moving  from 
one  to  an  adjacent  category)  occurs  between  class  D  and  class  C. 
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Figure  4.   Sigma-y  relative  vs.  downwind  distance  with  solid  line 
representing  the  recommended  parameterization  of  table  3.   Symbol 
A  represents  one  data  point,  B  represents  two  data  points,  etc. 
Same  data  as  figure  1,  except  long  range  data  is  included.   Data 
is  a)   Pasquill-Gif f ord  stability  class  B,  GMGO  data,  and   b)  class 
C,  GMGO  data. 
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Figure  5.   Same  as  figure  4,  except   a)   Pasquill-Gif f ord  stability 
class  D,  GMGO  data,  and   b)   class  E,  GMGO  data. 
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Figure   6.      Same   as   figure   4,    except      a)      Pasquill-Gif f ord    stability 
class   D,    NPS   data,    and      b)      class   E,    NPS   data. 
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DATA 

PASQUILL-GIFFORD 

CLASS 

SLOPE 

//POINTS 

GMGO 

B 

.0^167 

132 

GMGO 

C 

.0357 

no 

GMGO 

D 

.0227 

603 

GMGO 

E 

.0178 

206 

NPS 

D 

.0178 

215 

NPS 

E 

.0228 

72 

GMGO 

B  +  C 

** 

.0410 

272 

GMGO  + 

NPS 

D+E 

** 

.021  1 

1096 

**Last  2  rows  are  suggested  values. 

Table  3.   Relative  diffusion  parameterization  for  overwater, 
surface,  medium  range  (1-12km)  releases.   Slopes  are  least 
squared  linear  regression  to  data.   Sigma-y  relative  = 
( Slope )x( downwind  distance). 

Pasquill-Gif f ord      Sigma  Wind  Direction  (deg)       #  Hours 
Stability  Class       1  min  ave 1  0  min  ave       Of  Data 

B  3-96  6.41  7 

C  2.95  4.59  11 

D  1  .82  2.92  114 

E  1.32  2.61  11 

Table  4.   Sigma  theta  as  a  function  of  P-G  stability  class. 
Sampling  rate  is  1  Hz.  Data  is  from  NPS  experiments. 


This  two-class  feature  can  also  qualitively  seen  in  figure  7 
where  turbulence  intensity  as  a  function  of  stability  (defined 
using  the  Monin-Obuknov  length  [L])  was  plotted  by  Schacher  et 
al  (1982).   These  data  show  a  step  increase  in  sigma  theta 
(turbulence  intensity)  for  the  1  minute  averages  at  roughly  L  = 
-40m.   1  hour  averages  are  also  plotted  for  a  comparison,  even 
though  1  minute  averages  are  much  more  applicable  to  relative 
diffusion  in  these  experiments.   Note  that  the  1  hour  values  do 
not  show  this  effect  and  are  4  to  5  times  larger  than  1  minute 
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values.   This  suggests  that  our  two-class  approach  may  break  down 
as  the  cloud  travel  time  increases.   We  therefore  do  not  recommend 
interpolating  our  results  to  distances  greater  than  10  km. 
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(olej)  \z   - 


ZO/L 


Figure  7.  Wind  direction  sigma  (sigma  theta)  vs.  z/L 
(z  is  the  height,  10m,  and  L  is  the  Monin- 
Obukhov  length).  X's  are  1  min  averages, 
and  O's  are  1  hr  averages.  Data  is  from 
NPS.  (see  Schacher,  et  al. ,  1932) 
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4.   MEANDER  (SINGLE  PARTICLE)  PARAMETERIZATION 

For  purposes  of  describing  meander,  information  on  the  motion 
of  the  plume  is  available  only  from  the  NPS  data  set.   Meander  was 
measured  by  calculating  the  distance  between  the  center  of  mass 
for  a  given  profile  and  the  axis  defined  by  the  mean  wind  vector. 
The  method  used  to  calculate  the  mean  wind  vector  was  crucial. 
After  experimenting  with  various  techniques,  we  decided  to  let 
the  average  travel  time  of  the  plume  define  the  proper  averaging 
period  for  the  vector;  this  travel  time  was  roughly  one-half  hour. 
Since  the  wind  was  measured  only  at  the  release  point,  tracer 
profiles  were  correlated  with  average  wind  vectors  in  half-hour 
bins  offset  by  15  minutes  to  account  for  plume  transport. 

The  NPS  experiment  was  designed  so  that  the  plume  was  sampled 
at  nearly  discrete  ranges  from  the  release  point.   This  sampling 
procedure  allowed  us  to  segregate  data  into  range  bins.   The 
collection  of  off  axis  center  of  mass  positions  for  a  given  range 
bin  represents  a  probability  distribution  for  the  meander  of  a 
plume  or  puff.   Each  experimental  day  produced  several  distribu- 
tions at  various  downwind  distances  from  the  source.   The  distri- 
butions were  then  analyzed  for  normality,  and  the  standard  devia- 
tions (which  represent  sigma-y  due  to  meander)  were  determined. 

The  data  were  examined  for  dependence  on  windspeed, 
stability  class,  inversion  height,  experimental  day,  and  time  of 
day.   The  day  to  day  variability  was  so  large  that  it  tended  to 
swamp  other  dependencies.   Figures  8-21  show  day  by  day  plots  of 
the  center  of  mass  distributions  interlaced  with  the  respective 
wind  direction  time  histories  and  spectral  plots. 
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Figure  9.      a)     Wind  direction  tine  series,   and  b)     power  spectra 
for  tracer  release  depicted     in  figure  8.      Time  series  data  are 
2  minute  averages  of  wind  direction  sampled  at  1  Hz.      Spectral 
points  are  smoothed  with  a  5%  sliding  average.      S(n)    is  wind  dir- 
ection power  spectral  density  at  frequency  n.     D  is  the  total 
sampling  time.     <DIR>is  the  mean  wind  direction. 
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Figure  11.      a)     Wind  direction  tine  series,   and  b)   power  spectra 
for  tracer  release  depicted  in  figure  10.  / 
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Figure  13.  a)  Wind  direction  time  series,  and  b)  power  spectra 
for  tracer  release  depicted  in  figure  12. 
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Figure  15.   a)  Wind  direction  tine  series,  and  b)  powsr  spectra 
for  tracer  release  depicted  in  figure  14. 
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Figure  17.  a)  Wind  direction  time  series,  and  b)  power  spectra 
for  tracer  release  depicted  in  figure  16. 
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Figure  19.     a)     Wind  direction  time  series,   and  b)     pover  spectra 
for  tracer  release  depicted  in  figure  18. 
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Figure  21.  a)  Wind  direction  time  series,  and  b)  powsr  spectra 
for  tracer  release  depicted  in  figure  20. 
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Individual  data  points  are  marked  with  integer  values  of  the 
measurement  hour  for  a  given  experimental  day.   Note  that  the 
mean  position  of  the  cluster  often  moves  off  the  centerline,  but 
data  are  distributed  evenly  over  the  course  of  a  day.   Large 
variations  in  the  meander  "envelope"  can  be  seen  from  day  to  day, 
and  the  wind  direction  time  series  differ  radically.   The  power 
spectra  also  vary  greatly  from  day  to  day.   Spectral  parameters 
are  summarized  in  table  5.   Spectral  gaps  (defined  in  the  table  5 
caption)  exist  in  most  situations,  with  slopes  of  approximately 
-2  in  the  low  frequency  part  of  the  spectra. 
There  appears  to  be  no  clear  relationship  between  the  total 
energy  in  the  turbulent  and  low  frequency  parts  of  the  spectra. 
Since  meander  is  driven  by  the  low  frequency  energy,  this 
discourages  attempts  at  parameterizations  of  meander  based  on 
surface  layer  scaling. 

As  stated  above,  center  of  mass  distributions  were  grouped 
by  experimental  day  and  range  bin.   Variations  of  the  position  of 
the  measurement  aircraft  dictated  range  bins  to  be  2  km  long. 
Relevant  statistics  are  given  in  table  6.   Examination  of  the 
Shapiro-Wilk  statistics  (Shapiro  and  Wilk,  1965)  indicates  that 
the  crosswind  distributions  are  represented  fairly  well  by  normal 
distributions.   This  suggests  that  Gaussian  models  will  be  valid 
even  for  dispersion  caused  by  meander,  given  the  proper  sigma-y. 
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Date 
6/21 /82 
6/22/82 
6/24/82 
6/25/82 
6/27/82 
6/28/82 
6/29/82 


CI 

assicali 

P 

eak 

10 

min 

20 

30 

20 

5 

10 

9 

20  min 

50 

50 

25 

15 

15 

9 


Slopeg 

-6.5/3 

-7/3 

-6.5/3 

-6/3 

-5.5/3 

-5.5/3 


1  Classical  peak  is  the  inverse  of  the  frequency  associated  with 
the  "turbulent  energy  producing"  subrange  maximum  at  the  low 
frequency  end  of  the  -5/3  slope  (inertial  subrange)  region. 

2  Gap  is  the  inverse  of  the  frequency  where  a  minimum  occurs  on 
the  low  frequency  side  of  the  "turbulent  energy  producing" 
subrange . 

3  Slope  is  the  log  (power)  vs.  log  (frequency)  relationship  in 
only  the  "low  frequency"  part  of  the  spectrum.   Gap 
frequencies  are  not  included. 

Table  5.   Summary  of  wind  direction  spectra  for  NPS  data  depicted 
in  figures  9b,  11b,  13b,  15b,  17b,  19b,  and  21b. 
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DATE 

RANGE(km) 

N 

MEANwm^ 

STD    DEVpf 

rrn     SKEW? 

KURTOS IS u      W^ 

21Jun82 

2 

25 

-293.0 

149.7 

-0.2861 

-0.3701 

0.9705 

21 Jun82 

4 

26 

-356.5 

255.6 

+0.6759 

+0.0095 

0.9554 

22Jun82 

4 

^3 

-1733.7 

392.1 

-0.0101 

-0.6845 

0.9721 

24Jun82 

2 

13 

-80.5 

110.3 

+0.8840 

+0.6898 

0.9077 

24Jun82 

4 

34 

-678.0 

391  .5 

-0.4567 

+1 .0885 

0.9590 

25Jun82 

2 

14 

-201 .3 

50.4 

+0.7334 

-0.8242 

0.8912 

25Jun82 

4 

38 

-740.9 

166.3 

+0.7936 

+0.8523 

0.9474 

27Jun82 

2 

9 

-185.9 

103.3 

-0.0621 

-0.0441 

0.9567 

27Jun82 

4 

17 

-895.4 

198.8 

-0.1541 

-0.1224 

0.9674 

27Jun82 

6 

4 

-1516.6 

415.4 

+  1  .2242 

-0.7663 

0.8893 

28Jun82 

2 

6 

-417.3 

83.4 

+0.8019 

-0.1 838 

0.9396 

28Jun82 

4 

12 

-795.0 

260.9 

-0.1653 

+1 .4008 

0.9470 

28Jun82 

6 

6 

-970.6 

160.5 

-1  .3268 

+1 .7325 

0.8868 

28Jun82 

8 

5 

-1165.7 

467.0 

+2.1248 

+4.6043 

0.6790 

29Jun82 

2 

3 

-613.9 

477.9 

+0.9695 

0.961 3 

29Jun82 

4 

12 

-913.0 

1008.9 

+0.6763 

+0.0379 

0.9374 

29Jun82 

6 

9 

-789.9 

966.7 

+0.8435 

-1 .0370 

0.8362 

29Jun82 

8 

1 1 

-1 177.6 

1380.2 

+  1  .2122 

+0.7888 

0.8509 

1   MEAN: 


STD  DEV 


3   SKEWNESS 


Mean  position  of  plume  center  (as  measured  from  the 
mean  wind  centerline)  of  mass  distribution  for  a 
given  range  bin. 

Standard  deviation  of  plume  center  of  mass  for  a 
given  range  bin. 

Measures  the  tendency  of  deviations  to  be  larger  in 
one  direction  than  the  other.   Values  are  unbounded 
Positive  values  indicate  a  "tail"  of  values  exists 
in  the  values  larger  than  the  mean,  negative  values 
indicate  a  "tail"  in  smaller  values. 


KURTOSIS:   Measures  the  heaviness  of  the  distribution  "tails" 
(-2.0  <  Kurtosis  <  °°).   Relatively  large  values 
indicate  normality  assumption  should  be  questioned, 
but  this  statistic  is  unstable  for  a  small  number 
of  samples  which  is  characteristic  of  this 
experiment . 


W: 


Shapiro-Wilk  statistic  for  the  center  of  mass 
distribution.   Measures  the  normality  of  the 
distribution  (0.0  <  W  <  1.0).   Values  close  to 
indicate  more  normal  distributions. 


1  .0 


Table  6.   Meander  statistics  for  NPS  meander  data  .   N  is 
number  of  profiles  for  a  given  crosswind  center  of  mass 
distr ibut  ion . 


the 
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While  we  were  not  able  to  parameterize  meander  with  mean 
surface  layer  meteorological  variables,  direct  measurements  of 
lateral  turbulence  intensity  did  correlate  well  with  the  observed 
meander  standard  deviations.   Taylor  (1921)  predicts  that  when 
the  length  scale  of  turbulence  is  much  larger  than  the  scale  of 
the  plume  (sigma-y  meander  in  this  case)  the  following  simple 
formula  will  apply: 


ym 


i  x 
v 


(5) 


where  ovm  is  the  single-particle  lateral  standard  deviation,  or 


'ym 


sigma-y  meander, 


iv   is  the  lateral  turbulence  intensity, 

x    is  downwind  distance. 
This  is  often  referred  to  as  the  "near-field"  approximation  to 
single-particle  diffusion.   For  application  to  the  present 
analysis,  we  define  the  lateral  turbulence  intensity  as 


<v  > 

0 


(6) 


where  <v^>  is  the  variance  of  the  crosswind  component  measured 

over  the  travel  time  of  the  tracer  (approximately  1/2 
hour)  and  averaged  over  the  entire  day. 
U     is  the  mean  wind  speed 
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Figures  22-28  show  equation  (5)  plotted  for  each 
experimental  day  with  the  standard  deviations  of  table  5  (sigma-y 
meander)   depicted  as  data  points.   Also  shown  is  a  least-squared 
linear  fit  to  the  data. 

Note  the  good  agreement  between  this  simple  model  and  the 
observed  data  in  almost  every  case.   The  key  to  success  is 
matching  the  travel  time  of  the  plume  to  the  RMS  window  in  the 
turbulence  intensity  measurement.   Also,  since  the  measured 
variance  is  highly  non-stationary,  it  is  necessary  to  average 
this  quantity  over  a  full  day. 

The  agreement  between  the  data  and  the  near-field 
approximation  is  also  supported  by  our  spectral  results. 
The  spectra  continue  to  increase  in  power  through  the  entire 
spectral  window  (4-8  hr  "^ )  which  suggests  that  turbulent  motions 
operating  at  these  very  low  frequencies  have  very  long  time 
scales  as  equation  (5)  assumes. 
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Figure  22.   Sigma-y  meander  vs.  downwind  distance.   Error  bars 
are  the  standard  deviations  of  downwind  distance  positions  for 
a  given  range  bin.   Upper  number  is  Shapiro-Wilk  statistic  and 
lower  is  number  of  data  points.   Solid  line  is  best  fit  to  data 
Dashed  line  is  equation  5.   Data  is  NPS,  21  June  1982. 
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Figure  23.   Same  as  figure  22,  except  NPS  data,  22  June  1982. 
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Figure  24.   Same  as  figure  22,  except  NPS  data,  24  June  1982. 
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Figure  25.   Same  as  figure  22,  except  NPS  data,  25  June  1982. 
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Figure  26.   Same  as  figure  22,  except  NPS  data,  27  June  1982. 
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Figure  27.   Same  as  figure  22,  except  NPS  data,  28  June  1982. 
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Figure  28.   Same  as  figure  22,  except  NPS  data,  29  June  1982 
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5.   CHEMICAL  WEAPONS  HAZARD  FORECAST  PROGRAM  MODIFICATIONS 

Part  of  the  current  research  work  is  to  modify  the  Chemical 
Weapons  Hazard  Forecast  Program  (CWHFP).   This  section  assumes 
the  reader  is  familiar  with  that  program  and  understands  the 
program's  design,  inputs,  and  operation.   EPG  modified  the   24 
May  1984  "test  and  evaluation"  version  of  this  program  supplied 
to  EPG  by  J.  Branum  (1984).   Later  versions  of  this  model  may 
exist,  and  the  following  modifications  can  be  easily  transplanted 
into  such  hybrids. 

The  proposed  modifications  specifically  make  use  of  the 
parameterizations  described  above  and  are  designed  so  as  to 
supply  the  user  with  a  better  understanding  of  the  character  of 
atmospheric  diffusion  processes,  and  hopefully,  more  information 
on  how  to  avoid  hazards.   The  previous  model  used  meteorological 
inputs  to  calculate  a  stability  category,  and  then  calculate  one 
hour  averaged  plume  dimensions  based  on  stability  dependent 
sigma-y  and  sigma-z  algorithms.   This  version  uses  the  same 
stability  scheme,  but  replaces  the  one-hour  sigma-y  functions 
with  the  relative  diffusion  parameterization  presented  above. 
Sigma-z  parameterizations  are  unchanged.   The  resultant 
isopleths  represent  either  "dosage"  values  from  an  instantaneous 
"burst"  release  or  a  concentration  "snapshot"  from  a  continuous 
plume. 
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In  order  to  present  the  meander  of  a  plume  (offaxis 
excursions  of  a  puff)  the  single-particle  parameterization 
developed  above  is  used  to  obtain  a  "meander  envelope."   This 
envelope  is  simply  the  2  sigma-y  position  of  the  plume  or  puff 
center  of  mass  distribution  representing- a  95%    probability  of 
impact.   This  envelope  is  superimposed  on  the  instantaneous 
isopleths  in  order  to  show  the  combined  impact  of  the  two 
diffusion  processes.   An  artificial  "ripple"  is  convolved  into 
the  instantaneous  solution  to  offer  the  user  a  visualization  of 
meander.   The  frequency  and  amplitude  of  this  ripple  has  no 
physical  basis,  and  is  only  installed  to  help  the  user  understand 
the  nature  of  the  meander  process. 

Figures  29-36  show  several  examples  of  model  output.   Labels 
display  lethal  dose  levels  on  the  instantaneous  isopleths  and  a 
footnote  explains  the  meander  envelope.   Instruction  manuals 
should  explain  the  arbitrary  nature  of  the  ripple  in  order  to 
avoid  false  conclusions  based  on  these  plots.   Note  the  change  in 
the  frequency  of  the  ripple  as  the  meander  envelope  shrinks  (see 
figures  29,  30,  31  or  32,  33,  3*0.   This  change  again  was 
arbitrary  ,  and  was  only  intended  to  suggest  the  possibility  that 
the  important  eddy  sizes  decrease  as  the  meander  decreases.   The 
amplitude  of  the  ripple,  again,  was  arbitrary  (selected  to  be  1 
sigma-y  meander)  and  suggests  that  the  instantaneous  position  of 
the  cloud  is  somewhat  random  and  rarely  reaches  the  meander 
envelope  bounds.   Also  note  the  change  in  the  instantaneous  cloud 
dimensions  with  a  change  in  stability  class  (figures  29,  30,  31, 
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LD0-GD 


'10 


RING; 


300  YARDS 


BEARINGS  IN  DEG  TRUE 


PLOT  FORMAT  #1 


TERRAIN    TYPE 

flEHN    MIND 

SIGMA  WIND  DIRECTION 

MEAN  HIP  TEMP 

MEAN  SEA-SURFACE  TEMP 
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Figure  29.   Chemical  Weapons  Hazard  Forecast  Program  output 
sample  —  unstable  atmosphere,  large  meander. 
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CHEMICAL  WEAPON  HAZARD  TORECRST  PROGRAM  -  9NP.P  XX  .  X 
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Figure  30.   Chemical  Weapons  Hazard  Forecast  Program  output 
sample  —  unstable  atmosphere,  moderate  meander. 
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CHEMICHL  WERPON  HRZRRD  FORECAST  PROGRAM  -  SNRP  XV.  X 
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£OUPCE  SITE  i  effect iw*> 

SOURCE  RATE 


OPEN 

-SEA 

ie 

►ITS 

FROM 

7 

DEG 

10 

DEG  C 

1? 

DEG  C 

70 

e 

METERS  Bf 

SED 

ON  ♦* 

B 

MODIFIED 

PASOUILL 

MK  1  1 

5-SIZE  BOMB 'MI 

SSILE 

POINT-BURST 

88 

KG 

INST 

iNTANEOUS 

275  DEG  TRUE 


♦*  DEFAULT  ** 


CONTOUR  LABEL 
■  DOSE-AGENT) 
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Figure  31.   Chemical  Weapons  Hazard  Forecast  Program  output 
sample  —  unstable  atmosphere,  small  meander. 
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CHEMICAL    IJERPON    HRZfiRD    FORECAST    PPOGRflM    -    SNAP    XX. Y 
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Figure  32.   Chemical  Weapons  Hazard  Forecast  Program  output 
sample  —  stable  atmosphere,  large  meander. 
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CHEMICAL  HERPON  HAZAPD  TOPECRST  PPOGRP.M  -  9NAP  XX  .  X 
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Figure  33.   Chemical  Weapons  Hazard  Forecast  Program  output 
sample  —  stable  atmosphere,  moderate  meander. 
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Figure  34.   Chemical  Weapons  Hazard  Forecast  Program  output 
sample  —  stable  atmosphere,  small  meander. 
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CHEMICAL  ICRPON  HRZAPD  T0RECR9T  PROGRRM  -  SNRP  XX. V 
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Figure  35.   Chemical  Weapons  Hazard  Forecast  Program  output 
sample  —  light  wind  conditions. 
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CHEMICAL  tJEAPON  HAZARD  rORETCRST  PROGRAM  -  &NRP  XX.  X 
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Figure  36.   Chemical  Weapons  Hazard  Forecast  Program  output 
sample  —  strong  wind  conditions. 
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vs.  32,  33,  34)  or  wind  speed  (figures  35  vs.  36)  as  would  be 
qualitatively  predicted  by  a  Gaussian  model.   This  model  requires 
one  additional  user  input  over  those  needed  in  the  previous 
version;  the  standard  deviation  of  the  wind  direction  used  in  the 
meander  calculations.   As  stated  in  the  previous  sections,  the 
averaging  time  must  be  the  cloud  travel  time  to  the  distance  of 
interest.   Also,  the  value  used  should  be  the  average  of  several 
such  standard  deviations  obtained  over  a  significant  fraction  of 
the  day.   The  sampling  time  is  not  critical  ,  1.0  -  0.1  Hz  would 
be  sufficient. 

Code  modifications  are  listed  in  figures  37-40.   Figure  37 
shows  both  the  turbulence  intensity  input  and  the  revised  sigma 
parameter  table  (relative  diffusion).   Figure  38  shows  where  the 
meander  envelope  is  calculated  while  figures  39-40  list  the 
plotting  routine. 
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a) 


420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 


Mean_wi  nd_sp_k  t =5 
Mean  wind  dir=275 


i 


!  KTS 

!  DEG  TRUE 

!  METERS 


M i  x_l ayer_ht  =0 
! 

!  *$$$$  NEW  1985  $$$$$$$$* 

Turb_i nt ens= 1 2*P I / 1 80     !  turbulence  intensity  (radians)  input 

I  $$ tt$SS SifSSSSSSSSfiSSSJ 

I 

DIM  Ml_ht_est_meth$C30] 

Ml  ht  est  meth$="**  DEFAULT  **" 


b) 


3790 
3800 
3810 
3820 
3830 
3840 
3850 
3860 
3870 
3880 
3890 
3900 
3910 
3920 
3930 
3940 
3950 
3960 


!  LOAD  ARRAY  CONTAINING  VALUES  OF  Ap- >Dp 

FOR  P0=1  TO  6   !  WHERE  P0  =  NUMERICAL  REFERENCE  TO  STABILITY  CLASSES 
!  <  1  =  "A"  ->  6="F"> 

READ  Matr i x_ap<P0) , Mat r i  x_bp(P0) , Mat r  i  x_cp(P0) , Matr  i  x_dp(P0 ) 
NEXT  P0 

$$**$$$  NEW  1985  ***$$$$*$ 


DATA  0,0,0,0  ! 

DATA  0.0410,1.00,0.32,0.75 
DATA  0.0410,1.00,0.32,0.70 
DATA  0.0211,1.00,0.16,0.65 
DATA  0.0211,1.00,0.10,0.62 
DATA  0,0,0,0  !  — 

$t*t$*$st$ tssssssssststsss 


-  CAT  "A"  (NOT  USED  OVER  THE  OCEAN) 
--  CAT  "B" 

--  CAT  "C"   sigmas  are  for  relative  diffusion 
--  CAT  "D"   in  lateral  direction  only 
--  CAT  "E" 

-  CAT  "F->"<NOT  DEFINED  OR  USED  BY  THIS  MODEL) 


Figure  37.   Abbreviated  listing  of  CWHFP  showing  (a)  turbulence 
intensity  input  and  (b)  new  array  of  sigma  parameters. 
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5710 
5726 
5736 
5740 
5750 
5760 
5770 
5780 
5790 
5800 
5810 
5820 
5830 
5848 
5850 
5860 
5870 
5880 
5890 
5900 
5910 
5920 
5930 
5940 
5950 
5960 
5970 
5980 
5990 

6000 
6010 

6820 
6030 
6040 
6050 
6060 
6870 
6080 
6090 
6100 
6110 
6120 
6130 
6140 
6150 
6160 
6170 
6180 
6190 
6200 


I 
Pt_vapor_nomix: !  CflLC  MODEL  FOR  POINT  SOURCES  OF  GASEOUS  AGENTS, 
AND  NO  REFLECTIVE  MIXING  LAYER  CAP 


*** 

DIM  Me 
DIM  Dat 
Poi  nt  s_ 

!  $$*$ 

i 


M_t  o_y 

! 

! 

FOR  K  = 
Part  i  al 
Max_ran 
Dat  a_se 

X_i  nc  r 

IF    X_i 

! 

FOR  J=l 
X=J*X_i 
Dat a_se 
S  i  g  m  a_y 
S  i  gma_z 
Part i al 
IF    Part 

Dat  a_se 

NEXT  J 
NEXT  K 
i 

!  **** 

RflD 

FOR  J  = 
X=J*X_ 
IF  X>M 
Meand_ 
Meand_ 
GOTO  J 


$*    NEW  1985  *$*$ 

and_x< 100) , Meand_y( 100) 

a_set_x(3, 100) , Data_set_y(3, 100)   !  100  points  vs.  20  pts 

per_set=100  !  in  old  version 

$$t$$t$t$$t$ttt$ 

ds=39.37'36    !  CONVERSION  FACTOR  USED  IN  EQUATIONS  BELOW,  DEFINED 
AGAIN  HERE  TO  PERMIT  TESTING  OF  PARTIAL  PROGRAM 

1    TO    Num_dat a_set s 

1 (K)=Source_si  ze (K) / <PI *60*Mean_wi nd_sp_m*Dose_val <K) ) 
ge_x(K')  =  <Part  i al_l (K )/ ( Ap(K)*Cp(K ) ) ) ~(  1 /(Bp(K )  +  Dp(K > )  ) 
t_y(K,0)=0 

ement  =Max_range_x<K  > /Poi  nt s_per_set 
ncrement=0  THEN  X_i nc rement =1   !  ADDED  TO  PREVENT  BLOWUP  ON  BAD  DATA 

TO  Poi  nt  s_per-_set 
ncrement 
t_x<K, J)=X 

x=Ap(K)*X^Bp(K)   !  ASSUMES  NO  INITIAL  CLOUD  SIZE  FOR  CONSERVATISM 
x=Cp<K)*X-Dp(K) 

_2(K)=2*L0G(Partial_l(K)/(Si  gma_yx*Si  gma_zx) ) 
ial_2<K)<0  THEN  Part i al_2(K )=0  !  PREVENT  POSSIBLE  ERROR  IN  LATER  SQR 

CALCULATIONS 
t_y<K, J)=Sigma_yx*SQR(Part  i  al_2<K) ) 


*  NEW  1985  %%%%%%%%% 

1  TO  Poi nts_per_set 

i  ncrement 

ax_range_x(Num_data_sets)*C0S(2*Turb_i ntens)  THEN  GOTO  Close 

y < J)=X*TAN(2*Turb_i ntens)      !  meander  envelope  is  based 

x(J)=X  !  only   on  turbulence  intensity  and 

100  !  selected  to  be  the  2  sigma  value  (95"/.) 


Close:  Meand_y( J )=SQR(Max_range_x(Num_dat a_set s)~2-X^2)  !  stop  envelope 

Meand_x(J)=X  !    at  maximum  range  and  close  with  partial  circle 
J100:  NEXT  J 

MAT  Meand_x=Meand_x*<M_to_yds) 

MAT  Meand_y=Meand_y*(M_to_yds) 

i 


Figure  38.   Abbreviated  listing  of  CWHFP  showing  calculation 
of  the  meander  envelope.   Points  per  data  set  may  be  adjusted 
for  speed  or  better  resolution. 
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81*0   i  *t$$*    NEW  1985  $*$*****$$$ 

8200   DIM  Scaled_x< 190 ), Seal ed_y( 100)    !  LIMITS  BASED  ON  "POI NTS_PER_SET" 

8216   DIM  Contour_labl_x(3>,Contour_labl_y<3>   !  LIMIT  BRSED  ON  NUM_DATA_SETS 

S220   ! 

8230   FOR  Point_num=0  TO  Poi nt s_per_set    '  PO I NTS_PEP_SET  preset  above 

8240   Sc  al  ed_x  (Poi  nt_num  )=Meand_x  (Poi  nt_nu«i)*Sc  al  e_f  ac  t  or  !  scale  meander 

82  50   8c al ed_y  Poi nt_num>=Meand_y (Poi nt_num>*Sc al e_f ac tor  !  envelope 

8260   NEXT  Point_num 

8270   GOSUB  Plot_meander 

8280   ' 

8290   '  **$$*$$$$*$$$$$$$$$$*$**$* 

8300   ' 

8310   FOR  Set_num=l  TO  Num_dat a_set s        !  NUM_DPiTFl_SETS  PRESET  TO  "3",  ABOVE 

8320   ' 

8330   FOR  Poi nt_num=0  TO  Po i nt s_per_set    !  POINTS_PER_SET  PRESET   ABOVE 

8340   Seal  £d_x(  Po  i  nt_num>=Dat  a_set_x ( Set_num ,  Poi  nt_num  >*Sc  al  e_f  ac  t  or 

8  350   Seal  ed_y  (Poi  nt_num>=Dat  a_set_y (Set_num,  Po  i  nt_num)*Scal  e_f  ac  tor 

8380   NEXT  Point_num 

8  370   GOSUB  Plot_contour 

8380   Cent  our_l abl_x( Set_num)=Sc  al ed_x<Poi  nt  s_per_set ) 

8390   ! 

8400    i  $***$$*$$  NEW  1985  ttttt S$t$$S $ S * S S 1 $**$$ t$tS St 

8410   RAD 

8420   R  =  Dat  a_set_x  ( Set  _num , Po l nt  s_per_set )  /Dat  a_set_x ( Num_dat  a_set  s ,  Po  i  nt  s_per_s 

et  ■  i  d i mens i onl ess  distance  from  source  C0,1] 

8430   Omega=40*PI/180/Turb_i nt ens  !  ARBITRARY  frequency  of  ripple  oscillation 

8440   A=Scal ed_x(Poi nts_per_set >*TAN(Turb_i ntens)  !  ARBITRARY  ripple  amplitude 

8450   Ri ppl e=A*SIN(Omega*PI*R>  !  offset  for  y  position 

8480   C ont our_l abl_y <Set_num )=Sc al ed_y (Poi nt s_per_set  >+Ri pp  1  e  !  add  offset  to  la 

bel  position 

8470   DEG 

8438   '  ****$$$$$$$$*$$$$$*$$*$$$«**$$***$$*$$$**$$$$$* 

8490   ! 

8500   NEXT  Set_num 

8510   GOTO  Ex i t_dat a_p1 ot 

8520   ' 

8530   '  SSSSt$t    NEW  1985  *$$$*$*$*$$$ $$***$*$$$«$$$$$$*   (NEW  SUBROUTINE) 

8540   ! 

8550  PI ot_meander :   !  SUB  to  plot  meander  envelope 

8560   LINE  TYPE  1 

8570  Max_overst r i ke=l  !  same  as  "old"  version 

8580   X_offset=0 

8590   i  TOP  HALF 

86O0  FOR  St r l kecount =0  TO  Max_overst r i ke 

8610  Y_off set =. 33*Stri kecount 

8620  X_off set =. 33*Str l kecount  !  strikecounts  thicken  line 

8630  MOVE  Scaled_x(0),Sca1ed_y(0) 

8640  FOR  K=0  TO  Po i nt s_per_set 

8650  DRAW  Seal ed_x(K)+X_of f set , Seal ed_y (K >+Y_of f set 

8660  NEXT  K 

8670   !  SCALED_X,Y  now  contain  meander  envelope 

8680  i  PLOT  "BOTTOM"  HALF 

8690  MOVE  Sealed_x(0),Scaled_y(0) 

8700  FOR  K>0  TO  Poi nt s_per_set 

8710  DRAW  Scaled_x(K>+X_offset , -(Seal ed_y(K >+Y_off set ) 

8720  NEXT  K 

8730   NEXT  Strikecount 

8740  RETURN 

8750   !  $*$$***$$$$$$$$$$$*$*$$$$$* 

Figure  39.   Abbreviated  listing  of  CWHFP  showing  meander  envelope 
plotting  routine  and  contour  label  locating  scheme. 
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8776        !     *t*$$t$    NEW     1985    *$$*$$$$$$     (REVISED    SUBROUTINE) 

8786    PI ot_contour:        !       plot    contours    of    the    scaled    data 

8790       RflD 

8800  LINE  TYPE  1 

8810   Omega=40*PI/180  Turb_intens  !  ARBITRARY  frequency  of  ripple  oscillation 

8820  ! 

8830  !  PLOT  EEC  INNING  WITH  "TOP"  HALP  OF  CONTOUR 

8840  !   (  OVERSTRIKE  STUFF  OMITTED  ) 

8850  MOVE  Seal ed_x<0> , Seal ed_y<0> 

8860  FOR  K=0  TO  Po i nt s_per_set 

8870   R  =  Dat  a_s e t _x  (.  S e t _n u m ,  K  > <-  D  at  a_s e  t  _x  (  N u m_d at  a_s e t  s ,  P o  i  n t  s_p e r  s e t  )  !  d  i  mens  i 

onless  distance  from  source 

S8S0   A=Scaled_x(K)*TAN(Turb_int ens >  !  ARBITRARY  amplitude  of  ripple  selected  to 

be  1  sigma 
8890   Ri ppl e=A*SIN<0mega*PI*R>  !  offset  in  y  direction 
89O0  DRAW  Seal ed_x<K) , Seal ed_y<K>+Ri ppl e  !  add  offset 
8910  NEXT  K 
8920 

PLOT  "BOTTOM"  HALF    ....  most  same 


as  "top"  half 


8930 

8940 

8950  MOVE  Seal ed_> <0) , Seal ed_y(8> 

8960  FOR  K=0  TO  Po i nt s_per_set 

8  970   R  =  D at  a_s  e  t  _x  C  S e  t  _n u m  ,  K  >  ■'  D  at  a_s  e  t  _x  <  N  u m_d at  a_s e  t  s ,  P o  i  n t  s_p e  r  s  e  t 

8980   A=Scal ed_x ( K >*TAN< Turb_i nt ens > 

8990   Ripple=A*SINCOmega*PI*R) 

9000  DRAW  Seal ed_x<K> , -(Seal ed_y<K)-Ri ppl e>  !  subtract  offset 

9010  NEXT  K 

9020  DEG 

9030  RETURN   !  END  OF  GOSUB  PL0T_C0NT0UR 

9040   !  f$***$«$$*$$$$$$*$$$$$$$$$$$$$*$$$$$ 


Figure  40.   Abbreviated  listing  of  CWHFP  showing  plotting  of 
hazard  contours  based  on  relative  diffusion  parameterization. 
"Wiggles"  are  a  result  of  a  SIN  wave  imposed  on  the  Gaussian 
plume  model  solution. 
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6.   SUMMARY  AND  FUTURE  WORK 

The  presented  work  successfully  parameterizes  relative 
lateral  diffusion  from  mean  meteorological  quantities,  verified 
via  two  independent  overwater  tracer  data  sets.   Single-particle 
diffusion  is  successfully  parameterized  only  through  direct 
measurements  of  the  lateral  turbulence  intensity.   Once  this 
quantity  is  known,  however,  the  parameterization  is  well  behaved. 

Future  work  should  attempt  to  correlate  this  lateral 
turbulence  with  large  scale  synoptic  features,  mesoscale 
phenomena,  radiosondes  profiles,  or  other  more  easily  measured 
quantities.   The  GMGO  researchers  are  presently  being  petitioned 
for  positioning  information  on  their  experiments  so  that  a  more 
thorough  verification  of  the  single-particle  parameterization  can 
proceed  in  the  future. 

The  model  is  presently  designed  for  only  medium  ranges,  but 
could  be  extended  to  greater  distances.   This  extension  would 
make  results  highly  dependent  on  inversion  height.   Some 
prognostic  or  diagnostic  estimate  of  this  quantity  could  be 
incorporated  into  the  model.   Longer  range  overwater  experiments 
should  be  researched  for  plume  paramterizations  at  these 
distances . 
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